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We formulate a time-optimal approach to adiabatic quantum computation (AQC). A corresponding natural 
Riemannian metric is also derived, through which AQC can be understood as the problem of finding a geodesic 
on the manifold of control parameters. This geometrization of AQC is demonstrated through two examples, 
where we show that it leads to improved performance of AQC, and sheds light on the roles of entanglement and 
curvature of the control manifold in algorithmic performance. 
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Introduction. — Quantum computation is most commonly 
formulated in the language of the "circuit model" [1]. The 
problem of finding optimal quantum circuits - which mini- 
mize the number of gates used - was recently addressed in an 
elegant differential-geometric framework, wherein the mini- 
mum number of elementary gates for construction of a gen- 
eral n-qubit unitary U can be found by traversing the geodesic 
connecting the identity 1 to U over the SU(2") manifold |2|]. 
This approach is appealing since it allows for the application 
of powerful tools and techniques from variational calculus and 
differential geometry [3] to quantum computation. Adiabatic 
quantum computation (AQC) |4|], on the other hand, is a very 
different approach which has recently attracted much attention 
due to its fundamental connection to quantum many -body sys- 
tems, in particular to quantum phase transitions (QPTs) [5]. 
The basic strategy of AQC is to solve computational problems 
based on adiabatic evolution. A quantum system is prepared 
in the ground state of an initial Hamiltonian H(Q) = Hj. The 
system is then driven to a final state, which - provided the 
evolution was adiabatic - is the ground state of a "problem 
Hamiltonian" H(T) = Hp. This final state corresponds to the 
solution of a hard problem, while a short time T corresponds 
to an efficient AQC strategy. Although AQC is equivalent in 
computational power to the circuit model ffj, and relies on 
one of the oldest theorems of quantum theory - the adiabatic 
theorem [7] - it is still relatively unexplored. Specifically, an 
optimal strategy for AQC, akin to what has been done for the 
circuit model yD, has not yet been formulated. 

In this work we reformulate AQC as a variational prob- 
lem and develop a time-optimal strategy - a "quantum adi- 
abatic brachistochrone" (QAB) - for quantum algorithms [8]. 
Specifically, we devise a variational time-optimal strategy for 
obtaining an interpolating Hamiltonian H (t) between Hj and 
Hp, which gives rise to the shortest time T while guarantee- 
ing that the actual final state (the solution to the correspond- 
ing Schrodinger equation), is close to the desired final ground 
state. We go further and show that the QAB can be recast in 
a natural differential-geometric framework. Specifically, we 
construct a Riemannian geometry, along with its correspond- 
ing metric, for the adiabatic evolution. We provide two exam- 
ples which illustrate the advantage of this optimal approach. 



Time-optimal AQC. — The adiabatic approximation, which 
underlies AQC, is often stated as follows [7]. Consider a sys- 
tem subjected to a time-dependent Hamiltonian H(t), with a 
nondegenerate ground state |3>o(X)) isolated by a nonvanish- 
ing gap A(t) from the first excited state |$i(t)). Let D{t) = 
\{$l(t)\d t H(t)\$o(t)}\. Prepare the system in |^(Q)) = 
|$o(0)) and let it evolve according to the Schrodinger equa- 
tion into the state \ip(T)), Then, provided the time variation 
of the Hamiltonian is sufficiently slow, or T is sufficiently 
large, in that max (6 [ T ] D(t)/ rnin tg [ 0i x] A 2 (i) <C e, the fi- 
delity F(T) = | ($ (T)\ip(T))\ between th e final g round state 
and the actual final state is high: F > \/l — e 2 - However, 
it is well known that the latter condition is not always accu- 
rate 09[], as recently verified experimentally 111 Oil . Rigorous 
versions of the adiabatic approximation [ 1 1] typically involve 
\\dtH\\ 2 / A 3 (with the norm being the maximum eigenvalue), 
or terms with different powers of 1 1 <9 t if || and A. As our ap- 
proach is to use the adiabatic condition as a heuristic for find- 
ing optimal trajectories, the exact form of the adiabatic con- 
dition is in fact not essential: we shall judge success by the 
tradeoff between fidelity F and evolution time T. As we show 
below, this pragmatic approach also allows us to find time- 
optimal and geometric formulations of AQC. 

The time-dependence of Hamiltonians usually comes from 
a set of control parameters x(i) = . . . , x M (t)) - 

e.g., electric or magnetic fields, laser beams, or any other ex- 
perimental "knob" - varying over a parameter manifold M., 
whence H = H[x(t)]. Varying the Hamiltonian for a given 
interval t e [t , ti] then translates geometrically into moving 
along a control curve (or path) x(i) in A4. We can reparam- 
eterize Ai via a dimensionless "natural parameter" s(t) [3], 
with s(0) = and s(T) = 1 (e.g., the normalized length), 
where v(t) = ds/dt > characterizes the speed by which 
we move along x[s(t)] E M. To make the adiabatic dynam- 
ics locally compatible with the geometric structure of A4, we 
modify the adiabatic condition into the following local (i.e., 
instantaneous) form lfl2ll : 

v(s)\\d s H(s)\\ 

— — <e tee [0,1]. (l) 

Hereafter the norm is the Hilbert-Schmidt (or Frobenius) 
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norm, defined as ||A||hs = \/Tr[^4tyl]. From the relation 
T = J Q ds/v(s), we define the adiabatic time -functional 

where x = d s x. Inspired by the local adiabatic condition £[), 
we choose the instantaneous "adiabatic speed" via the ansatz 

v. dd ( S ) = eA 2 ( S )/\\d s H(s)\\, (3) 

and hence - using Einstein summation - the Lagrangian 
£[x(s),x(s)] = ||i i a i iJ[x(s)]||/eA 2 [x(s)], where d t = 
d/dx l . This ansatz is sensible, in that adiabaticity is hindered 
when the gap closes, while it is favored when the variation of 
the Hamiltonian is slow. To simplify the analysis, from now 
on, we take the Lagrangian to be C! = C 2 ; this corresponds 
to a reparameterization which leaves the length of the curve 
solving the Euler-Lagrange (EL) equations invariant |[2|, |3J . 

T measures the time taken to traverse the curve x(s) from 
start to finish, subject to the local adiabatic condition ([TJ. Our 
goal is to minimize this time and thus obtain the time-optimal 
curve, or set of time-dependent controls. This optimal curve is 
the QAB. From variational calculus, the optimal path xqab (s) 
should satisfy <ST[x(s)]/<5x(s) = 0, which gives rise to the 
EL equations dd^C [x, x] /ds — d x C [x, x] . 

Some remarks regarding the QAB are in order, (i) The total 
real evolution time T (for a given fidelity F) is not necessar- 
ily the same as the time-functional T. The correct interpreta- 
tion is this: after finding the optimal path xqab(s) we solve 
the Schrodinger equation iv a dd s \ip) — H\ip) with H and v- d d 
computed along the optimal path, to find |-0[xq AB (£)]). The 
actual adiabatic error is then 5(T) = yl — Fqab(T) 2 , where 
Fq A b(T) = |<$o(T)|V[xqab(T)])|, for a given T. With this 
interpretation it is safe to take e = 1 in Eq. Q. Thus, we use 
T as a guiding principle to obtain optimal paths; the actual 
adiabatic time and error should be calculated independently, 
as per the prescription above, (ii) Equation ([TJ and our corre- 
sponding choice of v a d are by no means unique. They merely 
represent a convenient ansatz for our subsequent analysis, and 
it is quite possible that a better ansatz involving a different 
combination of \\d s H || and A exists, (iii) The presence of A 
in C implies that in order to apply our method, one needs to ei- 
ther find the gap exactly (which could be as hard as solving the 
problem itself), or restrict the interpolation to forms for which 
an explicit functional form for the gap can be obtained (e.g., 
in exactly- or almost-exactly-solvable models), or follow the 
scheme suggested to translate a quantum circuit model com- 
putation to an AQC (by which an interpolating Hamiltonian 
with an easily calculable gap is derived) [6], or estimate the 
gap by means of other methods (e.g., experimentally), (iv) 
The bottleneck of adiabatic algorithms is at the finite-size pre- 
cursor of QPTs where the gap becomes small (and vanishes 
in the thermodynamic limit) 15|] . The subtle point, however, is 
that due to the 1 1 x % di H \ | factor in the numerator of C, in princi- 
ple, there is the possibility that one can (at least partially) sup- 
press the vanishing gap and associated QPT [13]. The QAB 



will inherently seek to identify such criticality-suppressing 
control strategies (relative to the specific adiabatic condition 
we have adopted here), if they exist. Indeed, we shall see 
this in our examples, below, (v) In general, not every choice 
for the norm || • || yields an analytic adiabatic velocity. Our 
choice of the Hilbert-Schmidt norm is made to ensure ana- 
lyticity, and to simplify our calculations. We shall see be- 
low that this choice also enables a geometric treatment of the 
QAB problem, (vi) Modeling AQC usually necessitates a pa- 
rameterization of the Hamiltonian. In general, though, this 
parametrization is not unique. For simplicity, consider the 
parametrization H(x) = ^\ irVj, where {o~i}^L x are time- 
independent, noncommuting, linearly-independent Hermitian 
operators, chosen in accordance with the underlying structure 
of the optimization or physical problem in question. For ex- 
ample, in a multi-qubit system, {e^} could represent (two-) 
local interactions in the form of the tensor product of Pauli 
matrices. It can be shown that in these "interaction coordi- 
nates," the EL equations read: 

x k + T k ij (x)x i x j = 0, (4) 

where 

T k l3 = 2 (C l3 C kl m - 5 k i d j A - S k j d i A) /A, (5) 

dj = (C)ij = Tr^o-j], andC ,J = (CT 1 )^. (vii) Including 
various physical constraints into the variational description of 
AQC is natural through the Lagrange multiplier method. 

Geometrization of AQC. — Motivated by the existence of a 
differential-geometric description for circuit optimization H, 
and by the resemblance of Eq. (0]l to a geodesic equation 
03J, we reformulate QAB in a differential-geometric language. 
This in turn endows AQC with a Riemannian manifold struc- 
ture. An immediate advantage of this geometrization is that 
it equips AQC with the powerful techniques and tools of dif- 
ferential geometry, and shows that the problem of finding the 
QAB can be viewed as belonging to geometric control theory 
[14]. Furthermore, geometrization allows treating AQC and 
the circuit model in a universal geometric setting |2], which 
may suggest a natural alternative to Refs. B6[] for proving the 
equivalence of AQC and the circuit model. 

The transition from the EL equations to a geodesic equation 
is possible when one can write £'[x, x] = gij(x)x l x : ' , where 
g (with matrix elements gif) is a differentiable and invertible 
matrix [3]. From the definition of C for the QAB, we then 
obtain 

gy(x) = Tr^x)^ (x)]/A 4 (x), (6) 

which is the sought-after metric tensor. In interaction coordi- 
nates, for example, g(x) = C/A 4 (x), while it is represented 
differently in other coordinates |3J|. 

In this framework, then, the QAB is equivalent to the 
geodesic over (.M,g), and Eq. © is the geodesic equation, 
in which r fc y = ^g ll (djgki + dtgji — 9;3jfc) are the connec- 
tion coefficients [where g l i = (g _1 )y]. Since g^j oc A -4 



3 




FIG. 1: (color online) Left: Final-time error S(T) for the RC and 2-d geodesic paths for the Grover search problem, for n = 6 qubits. 
Squares (cyan) indicate where the 2-d geodesic path outerperforms the RC path (5 gco dcsic < <5rc); circles (red) correspond to the opposite case. 
Oscillations are due to 5(T) itself being highly oscillatory, with an envelope oc 1/T. Middle: 5(T) for the 2-d and 4-d geodesic paths, for 
n — 1. Cyan squares (red circles) indicate where the 4-d (2-d) geodesic path results in a smaller error. Right: Component -R1212 {x 1 ,x 2 ) of the 
curvature tensor for n = 3. The curves on the curvature surface show the critical line (vanishing gap asm 00), the RC interpolation, and the 
2-d geodesic (QAB). -R1212 is the only independent component of the curvature tensor in this case [Rijki = 2iii2i2 (ffifeflji — 9u9jk)/ det g]. 



(if the numerator does not contribute a power of A) we find 
r ~ g~ 1 dg ~ A^ 1 9A. Using standard expressions yj, 
one can calculate the Riemann curvature tensor R from the 
connection coefficients and the metric tensor, yielding R ~ 
d 2 g + gT 2 ~ A~ 6 . 

Examples. — We consider the following two-dimensional 
(2-d) interpolating Hamiltonian: 

H (^(s),* 2 ^)) = x 1 (s)P^ + x 2 (s)P^, (7) 

where = 1 — |a)(a| for the normalized vector |a) e H 
with dim(7Y) = N (similarly for f^), ccq = (a|b) is 
a known function of N alone, and a; 1 (0) = x 2 (l) = 1; 

= x 2 (0) = 0. We can always find |a^) such that 
|b) = «o|a) +ai|a J -), where (ala- 1 -) = 0, and ot\ = (a^b). 
Completing {|a), ja -1 -)} to a basis for Ti we can easily diag- 
onalize the Hamiltonian (jTj, and find that the gap between 
the ground state and the first excited state is A(x 1 ,x 2 ) — 
VO 1 ) 2 + O 2 ) 2 + 2(2|a | 2 - \)x x x 2 . While the general 2- 
d problem requires a numerical solution, we find that if we im- 
pose a one-dimensional (1-d) constraint, i.e., x(s) = x 2 (s) = 
1 — x 1 (s), then an analytical solution is possible: 

£qab(s) = i a ° tan [(1 - 2s) arccos \a \] . (8) 

1 2^/1 -\olqY 

We now consider two illustrative problems which are special 
cases of the Hamiltonian (Q. 

Quantum search. — As a first illustration we revisit Graver's 
unstructured search problem, which involves finding a marked 
object among N objects by repeated oracle queries [1, 
Il5n . Graver's quantum circuit model solution uses 0{-/N) 
queries, which is provably optimal, and a quadratic improve- 
ment over the best possible classical strategy. This problem 
was successfully recast in the AQC setting by Roland and Cerf 
(RC) fl2Tl . who considered the 1-d version of (0 with x(s) = 
x 2 (s) = 1 — x 1 (s), I a) = J2k=o |fc)/V^V (equal superposi- 
tion), |b) = \m), and the fixed index m £ {0, . . . , N — 1} 
being the "marked item". Thus a Q = 1/VN, with N = 2 n 
the dimension of the Hilbert space of n qubits. It turns out 
that the optimal 1-d solution ([H) coincides precisely with the 



solution found by RC, who proved its optimality (in the sense 
of 0(yN) scaling for a fixed error) without the use of vari- 
ational optimization. We now extend the analysis by con- 
sidering 2-d and 4-d parametrizations, which corresponds to 
finding optimal curves on 2-d and 4-d manifolds, respectively. 
The 2-d case is given by Eq. (|7]) and the discussion that follows 
it, with I a) and |b) as above. In the 4-d case, we first con- 
sider a general one-qubit Hamiltonian H (x 1 , x 2 , x 3 , x 4 ) — 
-^(x 1 1 +x 2 (T x + x 3 a y + x 4 a z ), and solve the corresponding 
geodesic (or QAB) differential equations. Next we recall that 
Graver's search is effectively a 2-d problem (in the { | a) , | m) } 
basis). This enables us to use the 4-d setting for finding a 
Groverian geodesic path, with the proper boundary conditions 
corresponding to |a) = y/ (N — 1)/N\Q) + 1/VN\1) and, 
for example, \m) — |1). However, note that the parametriza- 
tion of H (x 1 , x 2 , x 3 , x 4 ) above is not the most general 4-d 
parametrization when n > 1. 

The 1-d RC analysis employed the local adiabatic condition 
([T} to recover the optimal scaling T opt oc V^V, for N 1. 
This might suggest that there is no room for further improve- 
ment, but we recall that in the AQC setting the fidelity is 1 
only in the limit T — > 00. Thus we compare the error 6(T) 
for the RC interpolation to the error obtained from the opti- 
mal 2-d interpolation. The result for n = 6 qubits is shown in 
Fig.[T](left); results for other values of n are qualitatively sim- 
ilar, though the advantage of the optimal interpolation shrinks 
as n grows. The optimal 2-d interpolation results in a smaller 
error for most values of T, a tendency that increases as T 
grows. Conversely, for most values of the error 8 the 2-d QAB 
requires a smaller time T than the RC curve. The middle panel 
shows the further improvement resulting from the 4-d interpo- 
lation. These results provide a rather striking demonstration 
of the power of our formalism, as due to its highly optimized 
nature, the Grover example is one where hardly any improve- 
ment was to be expected. 

Figure[T](right) depicts the RC and 2-d optimal curves over 
the curvature i?i2i2 surface. Clearly, the optimal curve fol- 
lows a path of lower curvature. This is confirmed in Fig. |2] 
(left), for different values of n. Figure [2] (right) depicts the 
amount of bipartite entanglement along the RC and 2-d opti- 
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FIG. 2: (color online) Red diamonds or solid line (blue circles or 
dashed line) represent the 2-d QAB (RC) path. Left: Average cur- 
vature component -R1212 vs n. Inset: instantaneous curvature for 
n — 4. Results for other values of n are qualitatively similar. 
Right: Same for negativity N12 11711 - between qubits 1 and 2. In- 
set: instantaneous gap (lower left) and negativity (upper right) for 
n — 4. Here, Ave 7 [Y] = J X[x(s)]ds/ f^ds for given path 7; 
X e {N12, -R1212} and 7 G {RC, QAB}. Note that same values of 
s (local time) do not necessarily indicate the same real (proper) time 
t for different paths. 



Conclusion and outlook. — We have presented a time- 
optimal, differential-geometric framework for AQC, and dis- 
cussed its implications for the optimal design of adiabatic al- 
gorithms. The power of this new framework was illustrated 
via an example showing how the performance of an adiabatic 
algorithm can be improved by increasing the dimension of the 
control parameter space, and how geometrization sheds light 
on the role of entanglement and control manifold curvature 
in this enhanced performance. The method presented here is 
general and can in principle be used to optimize any adiabatic 
quantum algorithm for which the gap (or estimate thereof) is 
known. An important next step is to incorporate decoherence- 
mi ligation strategies 12111 . 
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mal paths. In spite of its improved performance, there is less 
entanglement along the 2-d optimal path than along the RC 
path, so that more entanglement does not always translate into 
higher algorithmic efficiency. We have verified (not shown) 
that the same picture emerges in terms of the entanglement 
entropy (or block entanglement) jTg] . The explanation for 
this lower entanglement along the QAB is that it has a larger 
instantaneous gap than the RC path [Fig. [2] (right)]. Indeed, 
it has been shown that for Graver's algorithm the entangle- 
ment entropy is small away from the finite-size precursor of 
the first-order QPT, but peaks near the critical point ]5|], and 
we have verified the same for the 2-d QAB. Finally, the reason 
that QAB follows a path with larger gap is that this is consis- 
tent with higher adiabaticity. By our previous scaling result 
R ~ A~ 6 , it is also consistent with lower curvature. 

Linear equations. — Solving linear equations of the type 
Ay = a, where A is a given (Hermitian) N x N matrix and 
a is a given vector, is a common problem. Recently a quan- 
tum algorithm was proposed in the circuit model that can ob- 
tain y for well-conditioned, sparse matrices in a time scaling 
as polylog(Af) 11811 . Here we consider the problem of find- 
ing y = A _1 a as one of oracular adiabatic state generation 
(m . To do so we let |b) = ^4.~ 1 |a) /| | A" 1 |a) || . The formu- 
lation given above for the Hamiltonian (O then applies. For 
concreteness we let |a) = (1, . . . , 1) T f\N and take A to be 
an N x A*" Toeplitz matrix whose first row and column are 
successive natural numbers, starting from 1. Toeplitz matri- 
ces have important applications in signal processing j20ll . and 
are not sparse. We then find ao = y/2/N and hence can 
deduce immediately - by analogy to the Graver case, where 
ao = i/l/iV - that the optimal 1-d interpolation will give rise 
to a run-time T scaling as 0{yHSf) for a fixed error. Moreover, 
2-d and 4-d interpolations will further improve the error at 
fixed run-time. It is interesting to note that the most efficient 
known classical algorithm for inverting an N x N Toeplitz 
matrix requires 0(N log 2 N) steps [20], though a direct com- 
parison is not possible due to our oracular setting. 
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